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Abstract 

We study the effect of new flavor changing SUSY phases arising in the squark mass 
matrix in semileptonic decays B -> X s £ + £~, and B -> K^£ + £~, {£ = e,/x). SUSY 
effects are parametrized using the mass insertion approximation formalism. Constraints 
on SUSY contributions coming from other processes (e.g. b — ► 57, B — > K*£ + £~) are 
taken into account. Chargino and gluino contributions to photon and Z-mediated 
decays are included and non-perturbative corrections are considered. We study the 
correlation between the CP (and forward-backward) asymmetries and the expected 
value of the inclusive branching ratio. Several possible scenarios are distinguished and 
discussed according to the mass insertion that is considered. 



1 Introduction 



Forthcoming 5-factories are assumed to provide a great amount of data on many decays of 
the B mesons which have not been observed so far A study of this decays both in the 
SM and in its extensions is compelling. 

The transitions B — ► X s £ + £~ decays have been deeply studied. The dominant pertur- 
bative SM contribution has been evaluated in ref. [[| and later two loop QCD corrections 
have been provided @-f|. Recently the two loop SM matching conditions has been worked 
out in ref. 0. cc resonances contributions to these results are included in the papers listed 
in ref. [f7|. Long distance corrections can have a different origin according to the value of 
the dilepton invariant mass one considers. 0(1 /ml) corrections have been first calculated in 
ref. |§ and recently corrected in refs. ||, |IIJ] . Near the peaks, non-perturbative contributions 
generated by cc resonances by means of resonance-exchange models have been provided in 
refs. [|, [II], [12|, |13| . Far from the resonance region, instead, ref. |L4[ (see also ref. fllS] ) 



estimate cc long-distance effects using a heavy quark expansion in inverse powers of the 
charm-quark mass (0(l/m^) corrections). 

An analysis of the SUSY contributions has been presented in refs. [|16|]-[[0| where the 



authors estimate the contribution of the Minimal Supersymmetric Standard Model (MSSM). 



In refs. |T7|, [18| the authors consider first a universal soft supersymmetry breaking sector at 



the Grand Unification scale (Constrained MSSM) and then partly relax this universality 
condition. In the latter case they find that there can be a substantial difference between 
the SM and the SUSY results in the Branching Ratios and in the forward-backward asym- 
metries. One of the reasons of this enhancement is that the Wilson coefficient C-j(Mw) 
(see section |2| for a precise definition) can change sign with respect to the SM in some re- 
gion of the parameter space while respecting constraints coming from b — *> 37. The recent 
measurements of b — > 57 J20] have narrowed the window of the possible values of Ci(Mw)- 



Hence, it is worthwhile considering B — > X s £ + £~ in a more general SUSY framework then 
just the Constrained MSSM. In reference [21j the possibility of new-physics effects coming 
from gluino-mediated FCNC is studied. Possible relevant SUSY contributions coming from 
neutral Higgs bosons loops in the large tan/3 regime are analyzed in ref. 5!|. The process 
B — > X s t + t~ in the flipped SU(5) SUSY-GUT model is considered in ref. [p^| . Effects of 
SUSY phases in models with heavy first and second generation sfermions have been recently 
discussed by the authors of ref. [P3[ . They study the effect of the fi and A t terms in the 
MSSM Lagrangian on the branching ratio of B — > X s £ + £~ and no flavor changing phase is 
assumed. Here we consider the impact of possible new phases in the flavor changing part of 
the squark mass matrix. 

A SUSY "model independent" analysis of these decays has been provided in ref. f25fl 
where SUSY effects are parametrized using the Mass Insertion Approximation (MIA) |26| . 
In this framework one chooses a basis for fermion and sfermion states in which all the 
couplings of these particles to neutral gauginos are flavor diagonal. Flavor changes in the 
squark sector are provided by the non-diagonality of the sfermion propagators. The pattern 
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of flavor change is then given by the ratios 

(Sf\ _ ^ m ii) 2 AB m 
\"ij)AB ~ 2 , (1) 

sq 

where (m{j)\ B are the off-diagonal elements of the f = u,d mass squared matrix that mixes 
flavor i, j for both left- and right-handed scalars (A, B =Left, Right) and M sq is the average 
squark mass (see e.g. [0]). The sfermion propagators are expanded in terms of the 5s 
and the contribution of the first term of this expansion is considered. The genuine SUSY 
contributions to the Wilson coefficients will be simply proportional to the various 5s and 
a keen analysis of the different Feynman diagrams involved will allow us to isolate the few 
insertions really relevant for a given process. In this way we see that only a small number of 
the new parameters is involved and a general SUSY analysis is made possible. The hypothesis 
regarding the smallness of the 5s and so the reliability of the approximation can then be 
checked a posteriori. Many of these 5s are strongly constrained by FCNC effects E7\, EBL 129 



or by vacuum stability arguments Nevertheless it may happen that such limits are 



not strong enough to prevent large contributions to some rare processes. An estimation 
of branching ratios and asymmetries was finally provided considering the limits on the 5s 
existing at the time. 



In ref. [31J it is noted that recent CLEO [ 32| (90% C. L.) upper limit on the branching 
ratio of 

BR (b° -> iT > + /i") < 4.0 x 10" 6 (2) 



gives limits on the MI stronger than the ones of ref. [^J and new constraints on Wilson 
coefficients are derived. 

So far in the SUSY-MIA analyses all 5s have been considered as real parameters as this 
is the relevant part in order to better contribute to branching ratios and FB-asymmetries. 
However eventual phases of these MI can modify other observables which measure the CP 
violation (CPV) in these decays. In this paper we want to investigate this possibility in 
the inclusive decay B — > X s £ + £~ and give an estimate also in the exclusive channels B — > 
K(*k + £~ for £ = e or /i. 

For the exclusive channels a comment is in order. While the SM estimate for the inclusive 
process B — > X s £ + £~ is quite well established, much more uncertain is the issue of the 
exclusive modes. Theoretical errors are still large for B — > K^£ + £~ due especially to the 
computation of hadronic matrix elements and/or form factors. A number of approaches 
has been used to determine them and the interested reader can look in ref. f3l| , |33f| and 
references therein. A generic analysis of possible new physics effect on these exclusive decays 



is presented in ref. ||34j| : however in this paper the authors do not make any hypothesis on the 
source of new physics. SUSY-MIA models can predict deviations from the SM much larger 
than the actual theoretical errors. From an experimental point of view the presence of new 
phases in the mass insertions is surely much easier to be observed in the exclusive channels 
than in the inclusive one. In order to give an estimate of CP violation in the exclusive 



channels we refer to the computation of the authors of ref. [31]. We find that the amount of 



CPV is practically the same on the inclusive and the exclusive transitions. The treatment of 



2 



the exclusive decay however allow us to take into account the constraint coming from Eq. 



and thus to update the results of ref. |25 



The measure of CPV gives us important informations on possible SUSY extension of 
the SM. Either the experimental measure will or will NOT reproduce within the theoretical 
error (say ~ 20% for the inclusive decay and ~ 30 — 40% for the exclusive ones) the SM 
expectation for the BR's and FBA's, the detection of CPV will represent an important 
test for supersymmetry. Moreover, in ref. pSl it was pointed out that all SUSY effects in 



B — > X s £ + £~ can be parametrized in terms of few Mi's, namely (S^^ll an d (<$23 ) lr- We 
have considered both the possibility that all the 5s contribute to the various observables and 
the case in which only the LL sector of the squark mass matrix is flavor changing. In all of 
these scenarios we have studied how the measures of BR's, FBA's and CPV are correlated. 
Given a measure of one observable it is so possible to deduce what are the expectations for 
the others. 

A very particular behavior is observed if only LL mass insertions are allowed. If these 
insertions are real, no big signal has to be expected. However if they acquire an imaginary 
part they can sizably enhance the values of the BR's. 

In all cases where an important increase of the BR's is provided the amount of CPV can 
be at best 2-3 times the SM M value although it can be of the opposite sign. 

Instead, a depression of the SM result can be explained only with the help of the LR 
insertions and higher values of CPV can be expected in this case. 

The paper is organized as follows. In order to estimate the CP violation a number of 
observables can be defined. We have considered CP asymmetries both in the branching ratio 
and in the FBA's and we have collected our basic definitions in section |^. In section 3 we 
summarize the existing constraints on mass insertions. Section 4 and 5 discuss the impact 
of SUSY on CP non- violating and CP violating observables respectively. Finally conclusions 
are provided. 



2 Operator basis and observables 

The effective Hamiltonian for the decay B — > X£ + £~ (X = X s , K, K*) in the SM and in the 
MSSM is given by 



eff 



AG F 
'72 



K* ts K tb 



.i=l 47F i=9 
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where 



(3) 



Ql = SLal^LaCLpYcLp, <?2 = S Lo rftJ) L pC L ^C La , 

Q3 = S La ^^b La QLpl^qLPi Qa = SLalfj.h/3 ^ QLfiK^QLai 

q=u,..,b q=u,..,b 

Q5 = S La ^^b La QR/3YQR/3, Qg = SLa1iib L /3 QRPl^QRa; 



Q9 = (s L j l J> L )l'fl, Qio = (s L <y ll b L )T'f'y B l , (4) 

(1 -y 5 ) 

K is the CKM-matrix and qL(R) = 1- We neglect the small contribution proportional 

to K* s K u b. This is the only SM source of CPV in these decays. Our approximation is thus 
equivalent to say that SUSY is the only origin of CP violation. This Hamiltonian is known 
at next-to-leading order both in the SM |4], [J and in the MSSM [118] , |i~9|j . The most general 
low-energy SUSY Hamiltonian also contains the operators 

Q' 9 = (wMl-fl, 

Qio = (sRjMl-fysl. (5) 

However the contribution of these operators is negligible and so they are not considered in 
the final discussion of physical quantities |3TJ] . SUSY contributions to other operators 



are negligible because they influence our observables at an higher perturbative order. 

The definitions of the differential branching ratio and of the forward-backward asym- 
metries are (s = (pe+ + p^-) 2 /m 2 for the inclusive process and s = (pe+ + pi-) 2 /M^ in the 
exclusive ones, 9 is the angle between the positively charged lepton and the B flight direction 
in the rest frame of the dilepton system) 

T(s) = , (6) 

ri d 2 T(B -> Xtt-) 

I dcosy — bgnfcost/j 

A ( \ J-i dcost'ds (r,\ 

Afb[s) = 71 d 2 r(B -> xin~) ' (7) 



A' 



dcos^ 
l-i dcosy ds 

/•1 d 2 T(B - X£ + r) 

/ dcoso' — bgnfcosi 

, , _ J-i d cos ds 

FB[) ~ a f, d 2 Y{B^X s tn 

dcoso 1 / ds 



1 J d cos 9 ds 

where X = X s , K, K*. Note that for the inclusive case it is usually considered the nor- 
malized ratio R(s) = T(s)/T(B — > X c ei>); the explicit formulae can be found in ref |]J. For 
what concerns the exclusive decays a discussion of the non-perturbative methods necessary 
to compute the form factors is beyond the scope of this paper. We have considered the 



approach of ref. |3T| of Light Cone Sum Rules and the theoretical errors considered in this 
reference. In this case SUSY contributions can be disentangled only if they are far beyond 
the form factors theoretical uncertainty. 

In order to compute CP-violating effects it is necessary to include the absorbitive parts 
of the matrix elements of the various operators. The only contributions come from Qg. In 
order to use only prescription independent quantities it is usually considered the following 
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effective coefficient 

Cf(s)=C e 9 ff (s) + Y ICS (s)- (9) 

Cg (s) includes both the mixing with of the operators Q\ — Qq and Qs and the pertur- 
bative estimation of the hadronic matrix elements. Its complete definition for the SM and 
MSSM can be found again in refs. [§], [|, F res (s) provides the Breit-Wigner ansatz for cc 
resonances 0, 

y» = V> E fc * ^^f rw ; do) 

with C (0) fc, = 0.875 (for i = 1,...,6) In the literature there are at least other three 

different parameterizations of the resonant cc contribution: the HQET-based approach ||, 
the KS-approach |0| and the LSW-approach 



Our final estimates of these observables include non-perturbative effects (0(1/ ml) 
HHChPT ||10|| , 0(l/m? c ) corrections) and we refer to refs. Jl|, [25| for their complete 
expressions. 

It is possible to construct several CPV observables using the branching ratios and the 
FBA's. We will consider the following. 

• The global branching ratio CP asymmetry, 

= rw-cPFW] 
/<i 3 {r( s ) + CF[r( s )]} 

• The local branching ratio CP asymmetry 

A?R(s) _ T(s)-CP[T(s)] 
Al [S) ~ T(s) + CP[T(s)] ■ 

• The FBA differences 

T>a,a> = A FB (s) (A' FB (s)) - CP [A FB (s) (A' FB (s))] . (13) 

In this case it is better not to normalize this difference with the sum of the two 
observables. Afb(s) and A FB (s) can have a different sign in different s regions and 
their sum can be nearly zero according to the value of the Mi's. 

All the effects coming from the mass insertion approximation can be included in formu- 
lae (fHI) writing the coefficients CV, C| ff (s), Ciq as 

(~i /~iSM I siDiag . (~iMI 

Cf(s) = (Cf(s)) SM + (Cf) Dia9 + (Cf) MI , 

Cio = + C™ a9 + C™ 1 (14) 

where all the contributions are evaluated at the Mb scale and the various C® ia9 summarize 
all the contributions coming from graphs including SUSY Higgs bosons and sparticles in 
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the limit in which we neglect all the mass insertion contributions (they would be the only 
SUSY diagrams if the scalar mass matrices were diagonalized by the same rotations as those 
needed by the fermions). The explicit expressions for C® ia9 can be found in ref. [18]. In the 
following we will use |^8[ [25] , 

f Cf a9 {M B ) = -0.18 

J C g*9( M B) = -0.35 • (15) 
{ C^ a9 (M B ) = -0.27 



Cf 11 have been estimated in ref. |?5] and within the same range of parameters (fi ~ —160, 
M g i ~ M sq ~ 250 GeV, M- t ~ 90 GeV, M„ ~ 50 GeV, M H ± = 100 GeV, tan/3 ~ 2) one gets 



Qs(M) 



16/23 



C 7 MI (M B ) = [ - 0.19(5 d 23 ) LL - 33.4(4^ - 1.75(<%) LL - 0.25(5^) LR ] (- a#( w>) 

C 9 M/ (M B ) = -1.2(5 2 U 3 ) LL + 0.69(<%) M - 0.51($ 3 ) LL 
C^Affl) = 1.75(^ 3 ) LL - 8.25(^ 3 ) Lfl . 

(16) 

In this equation M is the scale at which all gauge bosons, top quark and SUSY particles 
are integrated out at the same time. The masses of this particles range between say 80 and 
250 GeV so that we have chosen M at an intermediate scale of about the top quark mass. 
Corrections to this approximation are included in the errors of our estimates. 

The sign and the value of the coefficient CV has a great importance. In fact the integral 
of the BR is dominated by the \Cj\ 2 / s and C-jCg term for low values of s. In the SM the 
interference between O7 and Og is destructive and this behavior can be easily modified in 
the general class of models we are dealing with. 

The big coefficient of (5% 3 )lr in C^ 1 (M B ) implies that the final total coefficient C W (M B ) 
can have a different sign with respect to the SM estimate. As a consequence of this, the sign 
of F-B-asymmetries can be the opposite of the one calculated in the SM. 

Finally the values of the physical constants we use are reported in table [IJ. 



3 Constraints on mass insertions 

In order to give an estimate of CPV it is necessary to take into account all possible constraints 
on mass insertions. 



The most relevant Ss interested in the determination of the Wilson coefficients C7, Cg 

d,u 

LL,LR- 



and C w are (#23^ 



Vacuum stability arguments regarding the absence in the potential of color and charge 
breaking minima and of directions unbounded from below |3D|| give 



'2M2 + 2M- 2 m 
l(*S) M |<- V ^ ' *2gL. (17) 

For M sq < 300GeV this is not an effective constraint on the mass insertions. 



A constraint on (S^^ll can come from the possible measure of AM Bs . 
In fact the gluino-box contribution to AMg s p9| is proportional to (5: 



d \2 
23/ LL 



instance ref. p9"|). A possible experimental determination of AMg s , say 

AM Bs < 30 ps" 1 

would imply that 



23 ILL 



< 0.5 



fsee for 



;is) 



(19) 



for squark masses about 250 GeV. Moreover the LL up- and down-squark soft breaking 
mass matrices, at scales higher than the EW-breaking one, are related by a Cabibbo- 
Kobayashi-Maskawa rotation 



(m; 



d \2 
sq) LL 



K\Ml)\ L K 



so that the limit (|19| ) would be valid for the up sector too: 



LL 



< 0.5 . 



(20) 



(21) 



Strong constraints come from the measure of B — > X^. The branching ratio of 
this process depends almost completely on the Wilson coefficients CV and C' 7 which are 
proportional respectively to (823) lr or rl and (^23)^^- The most recent CLEO estimate 
of the branching ratio for B — > X s 7 is [f20 



B exp (B -> s 7) = (3.15 ± 0.35 ± 0.32 ± 0.26) • 10" 



(22) 



where the first error is statistical, the second is systematic and the third comes from 
the model dependence of the signal. The limits given at 95% C.L. are [20fl: 



2.0 10^ 4 < B exp (B -> s 7) < 4.5 10" 



(23) 



We can define a Cj (M B ) as 



Cf{M 1 



B, 



B exp (B -> s 7) 



KtsK tb 

K 



cb 



QaF 

7cg(z) 



(24) 

|) we 

(25) 



where F can be found for instance in ref. [36]. Considering the limits in Eq. 
find at 95% C.L. 

0.29 < \Cf(M B )\ < 0.41 . 

Actually \Cf(M B )\ 2 = \C 7 (M B )\ 2 + \C 7 (M B )\ 2 and the constraint given in Eq. @ 
should be shared between the two coefficients. However in order to get the maximum 
SUSY contribution, we observe that in physical observables C' 7 does not interfere with 
C 7 , the C 7 C 9 term is suppressed by a factor m s /m b with respect to the C7C9 one and 
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CjC'g is numerically negligible (in fact Cg is much smaller than Cg). For these reasons 
we choose to fill the constraint of Eq. ( f25"l) with Cy(Mb) alone. 

Thus, we can choose the total C? (Mb) anywhere inside the allowed region given in 
Eq. ( ^5|) still remaining consistent with the MIA. 

The limit we get for (S^lb, is of order 10~ 2 and this rules out Z-mediated gluino 
penguins contributions to Cg and C w . 

For what concerns (5% 3 )ll we fi n d that the constraint changes significantly according 
to the sign of C^ (Mb)- In this case it is important to consider both the positive and 
negative region as this delta can give a non negligible contribution to Cg and C\q. The 
limits depend on the choice of the parameters in the chargino sector; the numerical 
results given below are computed for M sq ~ 250 GeV, /i ~ —160 GeV, mi R ~ 90 GeV, 
My ~ 50 GeV, tan (3 ~ 2 (in ref. p5[] these are the conditions under which we find the 
best SUSY contributions). 



Constraints on mass insertions are provided by the updated CLEO [B2] upper limit on 
the exclusive decay B° — > K*fi + fi~ written in Eq. (^), while no new limits are provided 
by 

BR(B + -> K + l + r) < 5.2 x KT 6 , 



BR (fi * 



K*e + e-) < 1.0 x KT 5 . (26) 



In ref. [[!!]] constraints on Cg, C10 are derived in the hypothesis of real mass insertions 
while we now consider complex 5s. In order to perform our numerical analysis taking 
into account this new limit we use the form factors of ref. [3T|. Although a better un- 



derstanding of the strong dynamics is compelling, we think that a clear signal of SUSY 
should overcome the theoretical uncertainties related to the calculation of hadronic 
matrix elements. 

Choosing the input parameters given in [^T|] in such a way to obtain the loosest limits 
on the Wilson coefficients, the result we obtain for the integrated BR (B° — > K*fx + fi~) 
is 

BR(5° — > K*fi + n~) = 10~ 6 [l.52 + 0.0399 |C 10 | 2 + 0.0399 |C 9 | 2 + 1.676 |C 7 | 2 — 

0.365 Re[Ci ] + 0.294 Re[C 9 ] + 0.197 Re[C 7 ] + 
0.0488 Im[C 7 ] + 0.0168 Im[C 9 ] + 0.273 Re[C 9 C 7 *]] (27) 

where the coefficients CV,Cg and C\o are defined as follow 



C 7 = C? ia9 {M b ) + C^ I (M l 



b 



Cg = (C^^iM,) + (C e g ff ) MI (M b ) , 

C 10 = C? ia9 (M b ) + C£'(M b ) . (28) 



S 



4 CP non-violating observables 



In this section we present the numerical results obtained for the maximum enhancement and 
depression with respect to the SM expectation of the BR and of the FB— asymmetries for 
the inclusive decays B — > X s £ + £~ and for the exclusive channel B° — > (K, K*)£ + £~ both for 
£ = e and £ = /i. For the inclusive case we update the estimates given in ref. due to the 
new limits on the 5s coming from Eq. (0). 



The relevant formulae for the above quantities can be found in ref. |25| where non- 
perturbative (0(1/ ml) and 0(l/m 2 )) corrections are provided as well; we refer to J25| for a 
complete list of references. 

A useful way to study the dependence of these observables on the mass insertions is to 
observe that the integrated BR and non-normalized F-B asymmetry are polynomial in the 
Wilson coefficients. It is thus possible to write the following master formulae: 



(BR e ^) Xs = l(r 6 [(9.79, 6.53) + 0.168 |Ci | 2 + 0.168 |C 9 | 2 + (41.5, 13.0) \C 7 
1.54 Re[Cio] + 1.10 Re[C 9 ] - (22.0,2.86) Re[C 7 ] + 
0.172 Im[C 7 ] + 0.0541 Im[C 9 ] + 1.30 Re[C 9 C 7 *]] , 

(BR e ^) K = 10" 6 [0.519 + 0.0133 |Ci | 2 + 0.0133 |C 9 | 2 + 0.0373 \C 7 \ 2 - 
0.122 Re[Cio] + 0.113 Re[C 9 ] + 0.189 Re[C 7 ] + 
0.00783 Im[C 7 ] + 0.00468 Im[C 9 ] + 0.0444 Re[C 9 C 7 *]] , 



(BR, 



( A NN \ 
\. A FB)K* 



10" 6 [(2.25, 1.83) + 0.0481 |C 10 | 2 + 0.0481 |C 9 | 2 + (5.90, 2.13) |C 7 | 2 - 
0.441 Re[Ci ] + 0.351 Re[C 9 ] + (-2.31,0.203) Re[C 7 ] + 
0.0584 Im[C 7 ] + 0.0196 Im[C 9 ] + 0.338 Re[C 9 C 7 *]] , 
2.14 - 0.0542 Im[C 10 ] - 0.467 Re[C 10 ] - 1.98 Re[C 7 C* ] - 
0.246 Re[C g C; ] + 9.09 Re[C 7 ] + 1.13 Re[C g ] , 
1.21 - 0.0211 Im[do] - 0.263 Re[C 10 ] - 0.488 Re[C 7 C* ] - 
0.0882 Re[C 9 C* ] + 2.24 Re[C 7 ] + 0.404 Re[C 9 ] 



(29) 

(30) 

(31) 
(32) 
(33) 



where the suffixes X s , K and K* correspond to the decays B 
coefficient Cj are defined as in Eq. (f2"§|). 

Some remarks on the above equations are necessary. 

The integrated normalized global asymmetries are given by 

Ni 



(X s , K, K*)£+£~ and the 



/ A' FB (s) 



A 



NN 
FB 



BR, 



(e>M) 



where 



N 



x s 




K* 



K, 



0.104 



N. 



K,K* 



r -l 
1.54 ps 



2 n 7T 5 



cb 



KKtbl 



f(z)k(z) 



3.30 10' 



3.70 10' 



(34) 

(35) 
(36) 
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(the functions / and g can be found for instance in refs. |5|, 37], z = M^/Mf^ 



There is no difference in the integrated non-normalized asymmetry for electrons and 
muons within an accuracy of 0.1 %. The asymmetry for B — > K£ + £~ vanishes identically 

Taking into account the discussion of sect. || and the expressions fll6|) for the coefficients 
it is possible to scan the allowed parameter space looking for the extremal values of the 
observables we are interested in. We present the results of this analysis drawing scatter 
plots of the two integrated asymmetries as a function of the integrated BRs. In this way 
the correlation between the expectations of branching ratios and asymmetries is explicitly 
worked out. We report the plots concerning both the inclusive and exclusive decays. 

We consider two possible scenarios. 

A. In Figs. [l]-0 we have allowed all 5s to vary at the same time. The range of possible 
values of the asymmetries is quite broad and in order to compare the SUSY result with 
the SM expectation we have reported in table ^| the central SM values for the same 
observables. 

Even for values of the branching ratio around the SM expectation, a sensible SUSY 
contribution to the asymmetries is possible. An enhancement of the asymmetry is 
particularly favored when also a branching ratio increase is provided. In the MIA it is 
however possible to have a strongly depressed branching ratio with respect to the SM. 
In this case the possible values of the asymmetries can both be largely positive and 
largely negative and, in any case, very different from the SM expectations. Note then, 
the strong effect of the constraint from the exclusive channel B — > K*fi + fi~: it is in 
fact responsible for the sharp cuts which are present for high values of the BR's and 
of the asymmetries. 

B. In Figs. [3|-§] we have considered what can be expected if only the 5ll are different from 
zero. This more restricted scenario can occur in several models and is quite peculiar. 
In this case it is still possible a change of sign of Cj if <5ll's are real and interesting 
results can in principle occur. Moreover we note that if the 5ll& acquire an imaginary 
part, the BR for the inclusive decay can increase considerably with respect to the SM 
but a "depression" is hardly possible. In this case the behavior of the FBA's is greatly 
modified, in fact they are non more allowed to change sign. This is due to the absence 
of the big (#33) lr contribution to C^ 1 . Sensitive enhancements of the asymmetries 
are anyway possible. 



In tables [3|-[1 we present the results and the values of the Wilson coefficients for which 
the extremal values are obtained. In these tables Cf are the total SUSY contribution to 
the Wilson coefficients at the Mb scale. 



5 CP violating observables 

The study of the CP- violating asymmetries defined in sect, ^follows the same guide-lines of 
the last section. Now, due to the presence of the resonances, we are forced to perform the 
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analysis of the integrated CP asymmetries into s regions far from resonances. 

First we consider the low-s region (1 — 6 GeV 2 ), see fig. |5|. Here we have quite high values 
of the differential branching ratios but the value of the CP asymmetries are just few percents. 
An analysis presented in ref. compute a BR CP-asymmetry of about — 0.19to;;|g% for 
the SM in this region. In our computation we have not considered the SM contribution to 
CPV as we expect that SUSY can overcome it. This is indeed the case on a relevant part of 
the parameter space. 

Then we take into account the high-s region (14 - 23 GeV 2 ), see fig. g. However for 
this part a comment is in order. The branching ratio in this zone of the spectrum is small, 
plagued with resonances (see fig. |) and making any measure on it represents an experimental 
challenge. Also from the theoretical point of view the presence of resonances makes it 
difficult to perform a rigorous estimate of CP observables. While a better comprehension 
of the resonant behavior of the BR is necessary, we present here a simple estimate of SUSY 
effects one can expect according to our present knowledge. CP asymmetry in this part of 



the spectrum was considered in |39| and found to be potentially large. We have checked it 
happens also in our framework. 

Our results are plotted in figs. g-§. In the case of exclusive decays we find very similar 
results so that we have not reported the corresponding graphs. Let us comment on the two 
possibles scenarios of the previous section. 



A. In figs. we have allowed all 5s to vary at the same time. The extremal values 
obtained in this case are in tables We find that the values of the CP-asymmetries 
are much more constrained in the case of an enhancement of the BR. In the case 
of a "depression" of the SM result larger values of CP asymmetries are possible. In 
particular, from the first plot in fig. |6| we note that 

• the extremal values of / Aq R can be in absolute value up to 6 times the one 
predicted in the SM; 

• for values of the BR around the SM result, it is possible that the CP asymmetry 
be up to 2-3 times the one predicted in the SM. 

Thus, the measure of CPV in these decays would represent an important issue for 
SUSY signals. So far a SM computation of J Af R and / V AA is not reported to 
our knowledge in literature, as their measure is experimentally more difficult to be 
performed. However we think that similar considerations are still valid. 

B. In figs. only the contribution of the 5ll's is permitted. In figs. §-|9| for (^)ll = 
the CP asymmetries are zero and the value of the BR for the SM + Diagonal con- 
tributions is obtained. If an enhancement of the BR is realized the value of the CP 
asymmetry is about the one obtainable in the SM. In this case a more refined study 
of the problem is necessary in order to see the interference between the SM and the 
SUSY contributions. 



11 



6 Conclusions 



A detailed analysis of SUSY contributions to CPV in semileptonic rare B decays has been 
performed using the mass insertion approximation. We have estimated the amount of CPV 
that can be expected given the present limits on the Mi's. Several observable have been taken 
into account. We have studied the correlation between the observables and the integrated 
BR for the inclusive decays. We have discussed the possible scenarios that can be realized 
according to the particular insertion one considers. The largest deviation with respect to the 
SM are of course expected when LR mass insertions are present. In the case in which only 
LL insertions are taken into account a detectable enhancement of the BR can be expected 
if the insertion has got an imaginary part. However in this case CP asymmetries are of 
the order of the SM ones and a more involved discussion about the interference of the two 
contributions has to be performed. 
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m t 


173.8 GeV 


m b 


4.8 GeV 


m c 


1.4 GeV 


m s 


125 MeV 


M B 


5.27 GeV 


a s (m z ) 


0.119 


l/a el (m z ) 
sin 2 9w 


128.9 


0.2334 



Table 1: Central values of physical constants used in the phenomeno logical analysis 



BR( 


e,X s ) 


9.6 10" 6 


BR(fi,X s ) 


6.3 10" 6 


A FB 


(e,X s ) 


23.0 % 


A FB (n,X s ) 


23.0 % 


A' 

FB 


(e,X s ) 


7.1 % 


A' FB (v,X s ) 


11 % 


BR(e, K) 


5.7 10~ 7 


BR(fi, K) 


5.7 10~ 7 


BR{ 


e,K*) 


2.3 10~ 6 


BR(fi, K*) 


1.9 10~ 6 


A FB 


[e,K*) 


17.0 % 


A FB (fx,K*) 


17.0 % 


A'fb 


[e,K*) 


19.8 % 


A! FB {n,K*) 


24.4 % 



Table 2: Central SM values for physical observables 
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(jSUSY 


(jSUSY 


riSUSY 
°10 


10 6 BR(e) 


M 
m 


22.2 
4.19 


0.734 - 0.086/ 
0.0429 - 0.0052/ 


-1.17-0.182/ 
-1.72 - 0.863/ 


-0.326 + 6.69/ 
4.23- 1.05/ 


W 6 BR(n) 


M 
m 


17.4 
1.72 


0.734 - 0.086/ 
0.0429 - 0.0052/ 


-1.17-0.182/ 
-1.72 - 0.863/ 


-0.326 + 6.69/ 
4.23- 1.05/ 


A FB 


M 
m 


33.5 
-33.6 


0.613 + 0.254/ 
0.567- 0.305/ 


-1.37 + 0.881/ 
-1.98 - 0.722/ 


0.149 - 1.89/ 
7.65 - 1.52/ 


A' FB (e) 


M 
m 


23.7 
-19. 


0.612 + 0.0328/ 
0.561 + 0.165/ 


0.522 - 0.88/ 
-2.05 + 0.664/ 


-2.9 + 0.911/ 
7.74+ 1.65/ 


A' F M 


M 
m 


27.3 
-27. 


0.678 + 0.154/ 
0.561 + 0.165/ 


-1.15 + 0.805/ 
-2.05 + 0.664/ 


-0.514- 1.46/ 
7.74+ 1.65/ 



Table 3: Inclusive decays B — > X s £ + £~ . Maximum (M) and minimum (m) values of the BR 
and of the F-B asymmetries (scenario A). Cf USY are the total SUSY contributions to the 
Wilson coefficients at the M B scale. 









(JSUSY 


r<susY 


riSUSY 
°10 


10 e BR(e,K) 


M 
m 


1.35 
0.0979 


-0.0265 - 0.0448/ 
-0.0359 + 0.0051/ 


0.665 - 0.101/ 
-1.72 + 0.872/ 


-2.13-5.75/ 
4.67- 0.25/ 


10 6 BR(n,K) 


M 
m 


1.35 
0.0976 


-0.0265 - 0.0448/ 
-0.0359 + 0.0051/ 


0.665 - 0.101/ 
-1.72 + 0.872/ 


-2.13-5.75/ 
4.67- 0.25/ 


10 6 BR(e,K*) 


M 
m 


5.44 
0.742 


0.734 - 0.086/ 
0.0429 - 0.0052/ 


-1.17-0.182/ 
-1.72 - 0.863/ 


-0.326 + 6.69/ 
4.23- 1.05/ 


10 6 BR(fi,K*) 


M 
m 


4.84 
0.421 


0.657 - 0.079/ 
0.0493 + 0.112/ 


-0.572 - 0.218/ 
-1.7-0.51/ 


-3.51 - 1.59/ 
5.59 - 0.501/ 


A FB (K*) 


M 
m 


25.3 
-25.2 


0.678 + 0.154/ 
0.561 + 0.165/ 


-1.15 + 0.805/ 
-2.05 + 0.664/ 


-0.514- 1.46/ 
7.74+ 1.65/ 


A' FB (e)(K*) 


M 
m 


33.3 
-29.5 


0.612 + 0.0328/ 
0.561 + 0.165/ 


0.522 - 0.88/ 
-2.05 + 0.664/ 


-2.9 + 0.911/ 
7.74+ 1.65/ 




M 
m 


36. 
-35.6 


0.678 + 0.154/ 
0.561 + 0.165/ 


-1.15 + 0.805/ 
-2.05 + 0.664/ 


-0.514- 1.46/ 
7.74+ 1.65/ 



Table 4: Exclusive decays B — > (K, K*)£ + £ . Maximum (M) and minimum (m) values of 
the BR and of the F-B asymmetries (scenario A). 
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rsusY 


rsusY 


r<susY 
Wo 


in 6 BB(e) 


M 


19.3 


n 647 + 265/ 


n 57 j_ n 377/ 

V 7 . '.J l | V 7 . " 7 1 1 A. 


-5 8 - 386/ 




ill 


O.c/O 




_n 979 — n rm t 


_zl _L n rni &f 


10 6 B B(u) 

X \J 1 J ± L- 1 Lai J 


M 


14.6 


n _l n 265/ 


n 57 _l n 377/ 


_5 s _ n 386/ 




111 


6 59 


_n D1 81 — fl 00779 7~ 


_n 376 — 01 1 1 T 


_4 «3 _i_ n 01 1 


Ann 


M 


32.4 


645 _ n 235/ 

\J • V7 -Lt_7 V 7 . 7 7 -£ 


568 - 335/ 

V7 • \J V7 W V ' '.J.J 7 -£ 


_5 70 + n 342/ 




m 


18.9 


-0.0742 + 0.0316/ 


-0.456 + 0.0449/ 


-4.75 - 0.0459/ 


A' FB (e) 


M 


22.8 


0.605 - 0.0741/ 


0.511-0.105/ 


-5.74 + 0.108/ 




m 


2.76 


-0.0741 - 0.0173/ 


-0.455 - 0.0247/ 


-4.75 + 0.0252/ 


A' F M 


M 


26.8 


0.605 - 0.0741/ 


0.511-0.105/ 


-5.74 + 0.108/ 




m 


4.88 


-0.0741 - 0.0173/ 


-0.455 - 0.0247/ 


-4.75 + 0.0252/ 



Table 5: Inclusive decays B — > X s £ + £ . Maximum (M) and minimum (m) values of the BR 
and of the F-B asymmetries for the LL-only case (scenario B). 









(jSUSY 


(JSUSY 


rsusY 
Wo 


10 6 BR(e,K) 


M 


0.927 


0.65 + 0.253/ 


0.574 + 0.36/ 


-5.8 - 0.368/ 




m 


0.48 


-0.0741 - 0.0173/ 


-0.455 - 0.0247/ 


-4.75 + 0.0252/ 


l0 6 BB(n,K) 


M 


0.925 


0.65 + 0.253/ 


0.574 + 0.36/ 


-5.8 - 0.368/ 




m 


0.478 


-0.0741 - 0.0173/ 


-0.455 - 0.0247/ 


-4.75 + 0.0252/ 


10 6 BB(e,K*) 


M 


4.65 


0.647 + 0.265/ 


0.57 + 0.377/ 


-5.8 - 0.386/ 




m 


2.2 


0.0545- 0.0218/ 


-0.272- 0.031/ 


-4.93 + 0.0318/ 


10 6 BB(fi,K*) 


M 


4.01 


0.647 + 0.265/ 


0.57 + 0.377/ 


-5.8 - 0.386/ 




m 


1.76 


-0.0741 - 0.0173/ 


-0.455 - 0.0247/ 


-4.75 + 0.0252/ 


A FB {K*) 


M 


24.3 


0.605 - 0.0741/ 


0.511-0.105/ 


-5.74 + 0.108/ 




m 


13.9 


-0.0742 + 0.0316/ 


-0.456 + 0.0449/ 


-4.75 - 0.0459/ 


A' FB {e){K*) 


M 


32.8 


0.605 - 0.0741/ 


0.511-0.105/ 


-5.74 + 0.108/ 




m 


13.7 


-0.0742 + 0.0316/ 


-0.456 + 0.0449/ 


-4.75 - 0.0459/ 


A' F M(K*) 


M 


35.5 


0.605 - 0.0741/ 


0.511-0.105/ 


-5.74 + 0.108/ 




m 


18.6 


-0.0742 + 0.0316/ 


-0.456 + 0.0449/ 


-4.75 - 0.0459/ 



Table 6: Exclusive decays B — > (K, K*)£ + £ . Maximum (M) and minimum (m) values of 
the BR and of the F-B asymmetries for the LL-only case (scenario B). 
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riSUSY 
7 


riSUSY 
9 


rsusY 

1U 


AW 1 flow) 


M 


1.31 


248 + 347/ 


-1 78 + 652/ 


5 32 - 206/ 




m 


-1.54 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


AE R (hish) 


M 


43.2 


-0 0359 + 0051/ 


-1 72 + 872/ 


4 67 - 25/ 




m 


-44.1 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


Da> (low) 


M 


0.00585 


0.0439 -0.188/ 


-1.81 + 1.09/ 


5.14-4./ 




m 


-0.00619 


0.027 + 0.144/ 


-1.82 - 1.17/ 


5.22 + 3.33/ 


D A > (high) 


M 


0.392 


-0.0146 + 0.198/ 


-1.47- 0.439/ 


4.69 -2.72/ 




m 


-0.401 


-0.0514 + 0.00682/ 


-1.69 + 0.49/ 


4.97 + 2.19/ 


Af R (low) 


M 


0.292 


0.248 + 0.347/ 


-1.78 + 0.652/ 


5.32 - 0.206/ 




m 


-0.344 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


Af R (high) 


M 


4.24 


-0.0359 + 0.0051/ 


-1.72 + 0.872/ 


4.67-0.25/ 




m 


-3.81 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


Da (low) 


M 


0.00131 


-0.00781 + 0.0518/ 


-1.52 + 0.0931/ 


5.21 - 3.01/ 




m 


-0.00143 


0.027 + 0.144/ 


-1.82 - 1.17/ 


5.22 + 3.33/ 


D A (high) 


M 


0.116 


-0.0514 + 0.0821/ 


-1.64- 0.785/ 


4.94- 1.56/ 




m 


-0.123 


0.0199 + 0.0504/ 


-1.76 - 0.777/ 


4.87+ 1.68/ 



Table 7: Inclusive decays B — > X s £ + £ . Maximum (M) and minimum (m) values of the CP 
asymmetries (scenario A) in the regions of low s (1-6 GeV 2 ) and high s (14-23 GeV 2 ) . 
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(JSUSY 


(JSUSY 


°10 


A RR (K*, low) 


M 


1.24 


0.248 + 0.347/ 


-1.78 + 0.652/ 


5.32 - 0.206/ 




m 


-1.45 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


A RR (K*, high) 


M 


40.3 


-0.0359 + 0.0051/ 


-1.72 + 0.872/ 


4.67- 0.25/ 




m 


-40.7 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


D A > (K*, low) 


M 


0.00758 


0.027 + 0.144/ 


-1.82 - 1.17/ 


5.22 + 3.33/ 




m 


-0.0072 


0.0439 - 0.188/ 


-1.81 + 1.09/ 


5.14-4./ 


D A , (K*, high) 


M 


0.383 


-0.0514 + 0.00682/ 


-1.69 + 0.49/ 


4.97 + 2.19/ 




m 


-0.374 


-0.0146 + 0.198/ 


-1.47- 0.439/ 


4.69- 2.72/ 


A RR (K*, low) 


M 


0.225 


0.167 + 0.362/ 


-1.68 + 0.568/ 


5.53 + 0.539/ 




m 


-0.263 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


A B L R (K\ high) 


M 


8.3 


-0.0359 + 0.0051/ 


-1.72 + 0.872/ 


4.67- 0.25/ 




m 


-7.91 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


D A (K*, low) 


M 


0.00136 


0.027 + 0.144/ 


-1.82-1.17/ 


5.22 + 3.33/ 




m 


-0.00126 


-0.00781 + 0.0518/ 


-1.52 + 0.0931/ 


5.21-3.01/ 


D A (K*, high) 


M 


0.0369 


0.0199 + 0.0504/ 


-1.76 - 0.777/ 


4.87+ 1.68/ 




m 


-0.0356 


-0.0514 + 0.0821/ 


-1.64- 0.785/ 


4.94- 1.56/ 


(K, low) 


M 


0.986 


0.315 + 0.385/ 


-1.55 + 0.885/ 


4. + 0.441/ 




m 


-1.08 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


Ag R (X, high) 


M 


33. 


-0.0359 + 0.0051/ 


-1.72 + 0.872/ 


4.67- 0.25/ 




m 


-34.1 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


Af R (K, low) 


M 


0.177 


0.315 + 0.385/ 


-1.55 + 0.885/ 


4. + 0.441/ 




m 


-0.195 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 


Af R (K, high) 


M 


7.01 


-0.0359 + 0.0051/ 


-1.72 + 0.872/ 


4.67- 0.25/ 




m 


-6.95 


0.122 - 0.271/ 


-1.67- 0.791/ 


4.63 + 0.0593/ 



Table 8: Exclusive decays B — > (K, K*)£ + £ . Maximum (M) and minimum (m) values of 
the CP asymmetries (scenario A) in the regions of low s (1-6 GeV 2 ) and high s (14-23 
GeV 2 ) . 
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Figure 1: Inclusive decay B — > X s fi + fi . Scatter plots of the Integrated F-B asymmetries 
(Afb and versus the BR (in units of 1CT 6 ). Scenario A. 
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Figure 2: Exclusive decay B — > K*fx + fi . Scatter plots of the Integrated F-B asymmetries 
(Afb and A' FB ) versus the integrated BR (in units of 1CT 6 ). Scenario A. 
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ure 6: Inclusive decay B — > X s i + i . Scatter plots of Aq R and Da> versus the integrated 
(in units of 10~ 6 ) in the low-s region. Scenario A. 
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ure 7: Inclusive decay B — > X s £ + £~. Scatter plots of an d -^A' versus the integrated 
(in units of 1CT 6 ) in the high-s region. Scenario A. 
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Figure 8: Inclusive decay B — > X s fi + fj, . Scatter plots of Aq R versus the integrated BR (in 
units of 10 _6 ) in the low-s region. Scenario B. 
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Figure 9: Inclusive decay B — > X s fi + fi~ . Scatter plots of Aq R versus the integrated BR (in 
units of 10~ 6 ) in the high-s region Scenario B. 
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